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ABSTRACT 

We present an 880 /im Submillimeter Array (SMA) detection of the submillimeter galaxy SXDF 850.6. 
SXDF 850.6 is a bright source {Sssofim = 8 mJy) detected in the SCUBA Half Degree Extragalactic 
Survey (SHADES), and has multiple possible radio counterparts in its deep radio image obtained at 
the VLA. Our new SMA detection finds that the submm emission coincides with the brightest radio 
emission that is found ~8'' north of the coordinates determined from SCUBA. Despite the lack of 
detectable counterparts in deep UV/optical images, we find a source at the SMA position in near- 
infrared and longer wavelength images. We perform SED model fits to UV-optical-IR photometry 
{u, B^V^ R^i\ z\ 3.6 /im, 4.5 /im, 5.8 /im, and 8.0 /im) and to submm-radio photometry 

(850 /im, 880 /im, 1100 /im, and 21 cm) independently, and we find both are well described by 
starburst templates at a redshift of z :^ 2.2±0.3. The best-fit parameters from the UV-optical-IR 
SED fit are a redshift of z = 1.87lo.o75 ^ stellar mass of = 2.5lo.3 ^ ^^^^ extinction 
of Av = s.oVi Q mag, and an age of 720l^?«° Myr. The submm-radio SED fit provides a consistent 
redshift of z ^ 1.8-2.5, an IR luminosity of Ljr = (7-26) xlO^^ L©, and a star formation rate of 1300- 
4500 Mq yr-^ These results suggest that SXDF 850.6 is a mature system already having a massive 
amount of old stellar population constructed before its submm bright phase and is experiencing a 
dusty starburst, possibly induced by major mergers. 

Subject headings: galaxies: formation, galaxies: starburst, cosmology: observations, galaxies: high 
redshift, submillimeter 
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Millimeter /submillimeter surveys have revolu- 
tionized observational cosmology by uncovering a 
substantial new population of mm/submm-bright 
dusty s tarburst galaxies at high redshifts (SMGs 
/e.g.. ISmaiL lyison fc Blainl Il997l : iBarger et all Il99 
^"^ U 120021: 



iHug hes et al.l 11998 1: 'Scot t et al.l [20021; [C oppin e t al.l 
2006; Greve et al.l [2 004: La urent et all 120051: Scott et^ 
1 20081: iPerera et al.l |2008|; [Austermann et al. 2009f 
Ixamura et al.l [20091 and see also Blain et al. 2002 for 
a review). The energy source of mm/submm emission 
is primarily from intense star formation activity, with 
star formation rates (SFRs) of 100 Mq yr~^ to several 
1000 Mq yr~^, and possibly partially from an activ e 
galactic nucleus (AGN) (e.g., Alexan der et al.l |2005[ ). 
While spectroscopic observations of radio-identified 
SMG s find a median redshift of 2.2 ([Chapman et al.l 
|2005[ ). several SMGs have now been found up to and 
beyond z = 4 (e . g.. ICapak et aTl [20081 : ICoppin e t al.l 
[20091 : iDaddi et al.1 12009[ ). There is a hypothesis that 



SMGs are progenitors of present-day ni assive ellipticals 
(e.g., Lillv et al.l[T999l : ISmail et al.ll2004f ). however, little 
is known about their evolution process. 

While multi-wavelength analysis is essential in order 
to understand the nature of SMGs, the coarse angular 
resolution of single dish telescopes prevents a precise 
determination of the exact optical/NIR counterparts. 
One of the most successful ways to pinpoint the loca- 
tion of the submm emissi on is to obtain high reso l ution , 
deep radio imagin g fe.g. Ilvison et aIl ll998L l2000L 120021 : 
ISmail et al.l l2000l : iBarger et al. 2000\ Although this 
technique reveals robust radio counterparts of ^50%- 
80% of submm sources (e.g., Ilvison et all I2005L 120071 : 
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IWagg et "ani2QQ9f ), sometimes multiple radio counterpart 
candidates are found for a source. 

The most accurate means of achieving high precision 
astrometry on the submillimeter emission is clearly to 
observe with high angular resolution at the wavelength 
of the original detecti on. In this re spect, the Submil- 
limeter Array (SMA; iHo et al.l [20041) has p roved to be 
a powerful instrume nt ("lono et al."20 06allbl: IWang et al.l 
120071 : lYounger et all [2007. 2008; Ivison et al.ll2008[ ). 

Here we present the results from the SMA observa- 
tions toward an 8 mJy submm source, SXDF 850.6, de- 
tected in the SCUBA Half Degree Extragalactic Survey 
(SHADES). SHADES has observed a large area of the 
sky (720 arcmin^) with high sensitivity (la ~ 2 mJy) 
with the purpose of obtaining a statistically significant 
unbia s ed sample of s ubmillimeter sources (Mortier et alj 
120051: ICoppin et all 120061 1. The regions covered by 
SHADES are divided between two fields, the Lockman 
Hole and the ?>nh^Tn/ XMM- Newton Deep Field (SXDF). 
SXDF 850.6 is a source in the SX DF with multiple 
optical, IR, and rad io counterparts ( Ivison et al.l 120071 : 
IClements et al.ll2008f ). but no established submm source 
identification. The strongest radio emission has no con- 
firmed optical counterpart, but the two secondary radio 
peaks both have apparent optical associations. 

§ 2 outlines the observational and calibration details, 
and the results are presented in § 3. In § 4, multi- 
wavelength data are described. The results of SED model 
fitting using the photometry are presented in § 5. In § 6, 
we discuss the nature of SXDF 850.6. A summary is pre- 
sented in § 7. Throughout the paper, magnitudes are in 
the AB system, and we adopt a cosmology with Hq = 70 
km s-^ Mpc-\ VLm = 0.3, and Qa = 0.7. 

2. OBSERVATIONS AND DATA REDUCTION 

SXDF 850.6 was observed on September 21 and Oc- 
tober 7, 2004, and on October 9 and 14, 2005, using 
a compact configuration with 7-8 antennas of the 
SMA. The phase center was positioned at the submm 
source centroid which is a(J2000) = 02^ 17^ 29.80^ and 
5(J2000) = -05° 03' 26'^00. The unprojected baseline 
length ranged from 23 m to 139m. The SMA correlator 
was equipped with 2 GHz total bandwidth in each side- 
band, yielding a total of 4 GHz bandwidth for contin- 
uum observations. A continuum channel was generated 
by vector averaging all of the channels after calibration. 
The SIS receivers were tuned to 345 GHz for the USB, 
yielding 335 GHz for the LSB. Interferometric pointing 
was checked at the beginning of the track and the point- 
ing offsets were usually within ~ ±5'' (15% of the pri- 
mary beam) for all antennas. We used an integration 
time of 30 seconds. 

The r aw SMA data was calibrated using the MIR 
package (|Scoville et al.llT993l ). Passband calibration was 
done using bright QSOs and planets observed during the 
track. Antenna based phase calibration was done using 
J0238+166 (1.02 Jy), J0423-013 (1.67 Jy) and J0132- 
169 (0.84 Jy). The fiux levels of all sources were nor- 
malized using the the quasar fiux estimates derived from 
the primary fiux standard Uranus. Imaging was carried 
out in MIRIAD (Sault et al. 1995). Maximum sensitiv- 
ity was achieved by adopting natural weighting, which 
gave a final synthesized beam size of 2'.'32 x 2'.'19 (P.A. = 
79.1°) and an RMS noise of 1.2 mJy. Because the source 
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Fig. 1. — The SMA contours overlaid on the VLA 21 cm 
image. The contours represent —3, —2, 2, 3, 4, and 5cr (where 
Icr = 1.2 mJy). The synthesized beam (2(^32 x 2(^19) is shown 
in lower left. The dark circle represents the approximate size of 
the SCUBA/JCMT beam (FWHM - 14^0- The dark cross, white 
crosses, and red circles are the AzTEC position (Ikarashi et al. i n 
prep.), three radio counterpart candidates in I vison et al.l ((20071 ). 
and MIPS 24 /xm sources in the SWIRE catalog, respectively. 



is close to zero in dechnation, the uv coverage is under- 
sampled in the north-south direction and the resultant 
synthesized beam has multiple sidelobes at the 30% level 
along the north-south direction about 14" away from the 
synthesized beam. The final map was made after adding 
the data from all four days and two sidebands, and cor- 
rected for the attenuation by the 35^' primary beam of 
the SMA. The astrometric accuracy is likely dominated 
by the low S/N of the image, and we assess ~0V4. 

3. RESULTS 

The synthesized map is shown in Figure [TJ A source 
was detected at about 6 sigma significance 8^2 north of 
the SCUBA coordinates. The source appears slightly 
elongated in the northeast-southwest direction but this 
is likely caused by the low S/N, and the source is entirely 
unresolved with the SMA beam. Based on a point source 
fit to the visibilities, the derived fiux is 6.9 ± 1.2 mJy, 
and the coordinates are a(J2000) = 02^ 17^ 29.79^ and 
(5(J2000) = -05° 03' 18V65. The derived fiux is consis- 
tent with the SCUBA 850 /im fiux of 8.15 ± 2.2 mJy 
(|Coppin et al.l [2006i) within the uncertainties of both 
measurements. 

4. MULTI- WAVELENGTH DATA 

Multi- wavelength images around SXDF 850.6 are 
shown in Figure [21 While an obvious optical counterpart 
is not seen in the images from u to bands, the NIR 
images represented in JHK and longer wavelength im- 
ages show a presence of a counterpart. The photometry 
of the counterpart and upper limits at the SMA position 
are presented in Table [H 

4.1. Radio and Millimeter 

llvison et al.l (|2007[ ) identify three radio counterpart 
candidates for SXDF 850.6 in their deep VLA 21 cm 
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Fig. 2. — Multi- wavelength images of SXDF 850.6 with the SMA contours (2, 3, 4, and 5a). The size of each image is 10^^ x 10^^ and 
north is up. From left to right: rgb image of MOSAIC 11/u, SuprimeCam/S, and V] rgb image of SuprimeCam/i?, i' , and z'; rgb image 
of WFCAM/ J, if, and K; rgb image of IRAC/chl (blue), ch2 (green), ch3 and 4 (red); MIPS 24 /im; VLA 21 cm. 



map based on th e probability analysis method of 
iDownes et all (HqSGI ). All three radio sources have corre- 
sponding MIPS 24 /im sources (Clements et al.'"2008), as 
shown in Figure [TJ The SMA observation reveals that the 
correct radio counterpart is the most distant source from 
the SCUBA centroid, and the brightest in 24 /im emis- 
sion. This demonstrates the effectiveness of the SMA 
in identifying SMG counterparts, especially in situations 
with multiple counterpart candidates. 

The radio coordinates derived from reanalysis of 
the VLA image (Arumugam et al. in prep.) 
are a (J2000) = 02^ 17^ 29.755^ (±0.005^) and 
d (J2000) = -05° 03' 18V40 (±0V08), with peak and in- 
tegrated flux of 99.8 /iJy and 100.0 ± 10.6 /iJy respec- 
tively. The close similarity of the peak and integrated 
flux suggests that the radio emission is unresolved with 
the r!79 X ^.'51 beam (P.A. = 6.1°). 

SXDF 850.6 is detected at 1100 /im with deboosted 
flux of 3.9 ± 0.5 mJy (Ik arashi et al. in pr ep.) us- 
ing the AzTEC camera (jWilson et al.l I2008D on the 
Atacama Su bmillimeter Telescope Experiment (ASTE; 
lEzawa et all f2004. 2008). The AzTEC coordinates are 
consistent with the SMA position within a 2a error cir- 
cle. 

4.2. Mid-IR 

Spitzer/IRAC and MIPS data are taken from the 
SpUDS archive. The Scr detection limits are 0.58, 0.89, 
5.7, and 5.3 /iJy in IRAC bands (3.6, 4.5, 5.8, and 
8.0 /im), and 36 /iJy in MIPS 24 /im band. The pho- 
tometry was performed using the SExtractor package 
(|Bertin fc ArnoutsI Il996[ ). The photometry in IRAC 
bands is conducted with 3^8 diameter aperture, and the 
photometric zero-points of the images are determined us- 
ing colors of earl y-type stars betwe en i^T-band and IRAC 
bands following iLacv et "all (|2005f ). In order to derive 
total flux, the aperture photometry are corrected by as- 
suming IRAC PSFs by a factor of 1.36, 1.40, 1.65, and 
1.84, respectively. The photometry in the 24 jam band 
is conducted using 12'' diameter aperture. The aperture 



magnitude is corrected by a factor of 1.70. 

We find emission at the SMA position in all IRAC 
bands and MIPS 24 /im band. Note that the 24 /im flux 
in Table [T] should be used as an upper limit since sources 
around the SMA position are blended in the 24 /im image 
(Figure [21). 

4.3. JHK 

J, H, and K band images are obtained from UKIDSS 
(UKIRT InfraRed Deep Sky Surveys) Ultra Deep Survey 
Third Data Release (jLaurent et al.l l2005). The 3a detec- 
tion limits are 24.28, 24.08, and 24.18 mag, respectively. 

A faint source is detected at the SMA position in all 
bands. We perform photometry using the SExtractor. 
Aperture photometry with 1'.'8 diameter aperture in the 
J, and K bands are corrected by difference between 
the aperture and the total (MAG_AUTO) magnitudes in 
the i^-band. 

4.4. uBVRi'z' 

We present a ix-band image of MOSAIC II on CTIO 
4-m telescope (Fujishiro et al. in prep.), 5, V, i?, z', and 
z' images of SuprimeCam on Subaru telescope in Figure 
[21 5, V, i?, z', and z' band images are obtained from 
Subaru/XMM-Newto n Deep Survey database (SXDS; 
iFurusawa et al.|[2'008f ) . No galaxy is detected at the SMA 
position in these bands. The 3cr detection limits are 
26.00, 27.94, 27.68, 27.43, 27.45, and 26.43 mag, respec- 
tively. 

Near the SMA position, three galaxies are seen; ~2" 
south-east, ~0'.'8 south-east, and ^1" south-west to the 
SMA position. We derive photometric redshifts for the 
sources and find that they are likely to be at z < 1, and 
it is unlikely that they are related to the SMA source 
(see §[5l). 

5. SED FITTING 

We perform SED model fitting to the photometry data 
from UV to radio to estimate photometric redshifts and 
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Table 1. Fluxes of SXDF 850.6 


Band 


Flux 


Instrument 


Ref. 


u 


<0.0148 /iJy 


MOSAIC II 


1 


B 


<0.0248 /iJy 


SuprimeCam 


1 


V 


<0.0315 /iJy 


SuprimeCam 


1 


R 
i' 


<0.0397 /iJy 
<0.0389 /iJy 


SuprimeCam 
SuprimeCam 


1 
1 


z' 


<0.0996 /iJy 


SuprimeCam 


1 


J 


l-40io.2o MJy 


WFCAM 


1 


±1 


Q oi +0.23 T 

'^•^-'^-0.22 A^^y 


WPP A A/T 




I\ 


17-1-0 OQ < 1 T^r 

y.i/ it u.zo /iJy 


A^/^TT'P A i\/r 
VV r AlVl 




3.6 /im 


OO O -1- O Q . . T-cr 

zy.U it U.o /iJy 




j 


4.5 /im 


QO Q J- 1 1 , , Ttr 

oy.o it i.i /iJy 






5.8 /im 


54.3 ± 1.9 /iJy 


IRAC 




8.0 /im 


38.2 ± 1.4 /iJy 


IRAC 




24 /im 


630 ± 30 /iJy 


MIPS 




70 /im 


<30 mJy 


MIPS 


2 


160 /im 


<200 mJy 


MIPS 


2 


450 /im 


<81 mJy 


SCUBA 


3 


850 /im 


8.15 ±2.2 mJy 


SCUBA 


3 


880 /im 


6.9 ± 1.2 mJy 


SMA 


1 


1100 /im 


3.9 ±0.5 mJy 


AzTEC 


4 


21 cm 


100.0 ± 10.6 /iJy 


VLA 


5 


0.2-12 keV 


<6.3 X 10~^^ erg cm^^ g-^ 


XMM-Newton 


6 



References. — (1) This work; (2) Flux limits in the SWIRE catalog 
(Surace et al. in prep.); (3) Coppin et al. 2006; (4) Ikarashi et al. in 
prep.; (5) Arumugam et al. in prep. (6) FLIX: upper limit server 
for XMM-Newton data provided by the XMM-Newton Survey Science 
Centre 



Note. 



Limits are 3a. 



other physical properties (stehar mass, F-band extinc- 
tion, star formation timescale, and metahicity). The 
photometric data are divided into two wavelength ranges: 
(1) from UV to 8.0 /im, where stellar emission is domi- 
nant, (2) from submm to radio, where thermal dust and 
synchrotron emission is dominant. 



5.1. 



IR 



UV - Optical 

SED fits are performed to the photometry ranging 
from UV to 8.0 jam {u^B^V^R^i' ^ z' ^ J^H^K ^ 3.6 /im, 
4.5 /im, 5.8 /im, and 8.0 /im) using the Hyperz code 
(|Bolzonella et al.l l2000f ). We use a synthetic spec- 
tral library of GALAXEV (Bruzual & Chariot 2003) for 
SEP templates. We adopt the Calzetti e xtinct ion law 
(jCalzetti et al.l I2000D and the iChabrierl (|2003[ ) initial 
mass function with lower and upper cutofP mass of 0.1 
and 100 Mq. Two star formation histories are assumed: 
(i) exponentially decaying SFR with a timescale of r (i.e., 
SFR (X e"*/'^), (ii) constant SFR. We calculate val- 
ues for each SED template with free parameters of red- 
shift {z = 0-6), stellar mass (M^), T^-band attenuation 
{Ay = 0.0-6.0), star formation timescale (r = 0.01- 
5 Gyr), and metallicity [Z = 1.0,0.4,0.2 Zq). The 
best-fit results are obtained with the exponentially de- 
caying SFR, z = 1.87t^;J^, = 2.5tg-J x 10^^ M©, 

iv = 3.0t?;o mag, and age = 720t2?o° Myr (Table [2] 



-0.07' 



S.Otii mag, and age = 720 



+ 1880 
-210 



and Figure [3j). The results suggest that SXDF 850.6 
has a large population of old stars, with large amount of 
dust obscuring the stellar emission. The stellar mass 
is comparable to those of massive star- forming ^alax - 
ies at z - 2 (e.g., Shapley et al. 2005; Erb et al.' '2006a) 
and consistent with stellar masses rec ently derived for 
other SMGs (-lO^MO^^ M©) (e.g., iDve et al.l 120081 : 




z= 1.87, Ay =3.0, age = 720 Myr 
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Observing Wavelength [|j,m] 



Fig. 3. Best-fit SED obtained in UV-optical-IR SED fit 
along with photometry data. Downward arrows represent 3a upper 
limits. 



iMichalowski et al.|[2QQ9l ). 

5.2. Submm - Radio 

Photometry at 850 /im, 880 /im, 1100 /im, and 21 cm 
is used for submm-radio SED fits. The MIPS photom- 
etry is also used to check for consistency. We use the 
following starburst SED models: Arp 220 (|Silva et al.l 
Il998f ). th e average SED of 76 SMGs with spectroscopic 
redshifts fMichalowski et alji2009|), and 105 SED mod- 
els of Chary & Elbaz (2001). We perform minimized 
fits with free parameters of redshifts and flux scaling fac- 
tors. We find the best-fit redshifts of z = 2.11q^ for 

Arp 220 (x^ = 0.3), z = 1.8t^J for th e average SMG 
(X^ = 0.4), and z = 2.5t^ j for an SED of lCharv &: Elb^ 
(|200lf ) (x^ = 0.09). The errors are 99% confidence inter- 
va ls. Note t hat there are other templates in the library 
of iCharv fc Elbad (|200lf ) that fit the photometry data 
with similar x^ values in the redshift range of z ^ 1-3. 
The derived redshifts are consistent with the result in the 
UV-optical-IR SED fit in the previous section. The IR 
luminosities are estimated to be Ljr = (7-26) x 10^^ Lq 
from the intrinsic IR luminosities of SED templates mul- 
tiplied by scaling factors for the best-fit SEDs. The best- 
fit SED models are shown in Figure |4] with photometry 
data from UV to radio. All of the SEDs overestimate 
the UV-NIR photometry, suggesting that stellar emis- 
sion is heavily attenuated by dust. This is consistent 
with the large extinction {Ay = 3.0) derived from the 
UV-optic al-IR SED fit. Finally, w e note that the SMG 
sample in 'Michalowski et al.l (|2009f ) is biased against op- 
tically faint objects. It is clearly important to take into 
account a heavily obscured SMG li ke SX DF 850.6 (or 
e.g., GOODS 850-5 in Wang et al.l [20091 AzTECl in 
Tamura et al. in prep.) to properly understand the over- 
ah picture of SMGs. 

6. DISCUSSION 

6.1. AGN Contribution 

The rest-frame 1.6 jam bump is clearly seen at 5.8 jam 
(Figure [3]), suggesting that the NIR emission detected 
in the IRAC bands i s star-formation do minated (e.g., 
IWeedman erall l2006l : iFarrah et all I2008D . The flatter 
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Table 2. Best-fit Results in UV Optical IR SED Fit 



z 

(1) 


(Mo) 
(2) 


Ay 
(mag) 
(3) 


Age 
(Myr) 
(4) 


(Myr) 
(5) 


Z 
(6) 


(7) 


-1 07+0.15 




n+0.3 

•^•^-i.o 


79PI+1880 


20 


1 


13.4 



Note. — The errors are 68% confidence intervals. (1) photometric redshift; (2) stellar mass; (3) V-band attenuation; (4) age; (5) star formation 
timescale; (6) metallicity; (7) value 



X 



105 

104 
103 
102 

101 
100 
10-1 

10-2 































Arp220 
SMG 
CEOl 
1 


z = 2.1 

z=1.8 1 

z = 2.5 



(e.^.. iKovacs et ani2QQ6t [Pope ( 
"Mic habwski l2009h . 



10-1 10° 101 102 103 104 
Observing Wavelength [|j.m] 



10^ 106 



Fig. 4. — Best-fit SEDs for three different models obtained in 
submm-radio SED fit. The SEDs are the average SED of 76 SMGs 
(jMichalowski et al. 2009) at z = 1.8, Arp 220 (Silva et al. 1998) at 
z = 2.1, and a starburst template of Charv &; Elbaz ( 2001) at 
z = 2.5. Photometry data from UV to radio are overplotted. 



spectral slope in the IRAC emission compared to AGN 
populations f avours a star formation as a domin ant heat- 
ing source ( Yun et alJ[2QQ8l: llvison et alJ l2QQ4l. This is 
also supported by the fact that SXDF 850.6 appears 
close to the starburst model trac ks in the S h o/S4^ 
vs. S24/S8.0 color-color diagram of llvison et all ([2001 . 



These facts suggest that SXDF 850.6 favours star for- 
mation as a dominant heating source and the physical 
quantities derived in the SED fit are less affected by the 
AGN. 

The 3<j upper limit on the rest-frame X-ray luminosity 
derived from 0.2-12 keV flux is 1.3 x 10^^ erg s~^ (assum- 
ing z = 1.87 and an effective photon index of F = 1.8). 
Since this value is not properly corrected for hydrogen 
attenuation, we do not constrain the AGN contribution 
from the X-ray luminosity. 

6.2. Dust Mass and Molecular Gas Mass 

Assuming the observed 880 jam flux is dominated 
by thermal dust emission, the dust mass can be de- 
rived as Md = 5'obs^i/[(l + ^)^d(^rest)^(^rest, ^d)] (e.g., 
iHughes et al! 1997), where 5'obs is the observed flux den- 
sity, i^rest is the rest-frame frequency, /^d(^rest) is the 
dust mass absorption coefficient, Td is the dust tem- 
perature, and B{i'rest,Td) is the Planck blackbody func- 
tion. We assume that the absorption coeffi cient varies 
as K,d oc and /3 lies between 1 and 2 (e.g, iHildebrandl 
[l983). We adopt /^d(125 /im) = 2.64 ± 0.29 m^ k^-\ 
the average value of various studies (Dunne et al.ll20Q3l ). 
Td = 30-50 K, and P = 1.5, the typical values for SMGs 



^tal]l2006l : ICoppin et al.l 

The dust mass is esti- 
mated to be Md = (4-9) xlO^ M© for SXDF 850.6 
at z = 1.87. This is consistent with previous work 
on SMGs (e.g.. iKovacs et al.l 120061: ICoppin et all 120081 : 
iMichalowski et al.ll2009[ [ 

By adopting a gas-to-dust mass ratio of 54, wh ich is 
an average value for SMGs in Kovacs et al.l (j2006f ). the 
molecular gas mass is Mgas = (2-5) xlO^^ Mq. 

6.3. Star Formation Activity and Nature of SXDF 850.6 

The UV-IR SED exhibits quiescent star forming ac- 
tivity dominated by old stellar components. The cur- 
rent SFR estimated from the SED fit is approximately 
zero. This is because the large part of the stellar 
mass was formed at the early phase of star formation 
and SFR decreased with time following exponential de- 
cay (r = 20 Myr for the best-fit SED). On the other 
hand, the IR luminosity of the best-fit su bmm-radio 
SED provides SFRir = 1300-4500 M© yr-^ (jKennicuttl 
I1998D . The SFR derived from 1.4 GHz radio emission 
is 1400 Mq yr-^ following t he equation of [Belli (l2003f ). 
and 1100 M© yr"^ following lYun fc Carillil (|2002h (with 
a spectral index of —0.8 and a radio-FIR normalization 
factor of 1). 

To see whether the observational data allow the coexis- 
tence of the old stellar SED with quiescent star formation 
and the dusty starburst SED, we create synthetic SEDs 
at UV-IR wavelength using the GALAXEV library. The 
synthetic SEDs are composites of an SED with the same 
parameters as the best-fit SED which is dominated by 
old stellar population and starburst SEDs with differ- 
ent parameters. We find that the composite SEDs with 
plausible parameters for starburst SEDs (e.g., SFR ~ 
2000 Mq yr-^ Ay = 3.0, and age = 1 Myr) are con- 
sistent with the photometry data (Figure [5j). Note that 
combining the SEDs increases the total stellar mass by 
only ~1%. 

Combining these results allows us to infer the nature 
of SXDF 850.6: it is a mature system with a large frac- 
tion of old stellar components, and currently experienc- 
ing a vigorous dusty starburst. The coexistence of old 
sta rs and a cur r ent st arburst in an SMG is suggested 
by "Wang et al." (2009) from a detailed SED analysis. 
SXDF 850.6 has enough molecular gas mass as estimated 
in § 16.21 to maintain intense star formation. Such sig- 
nificant s tar formation is likely caused by major merg- 
ers fe.g.. iGreve et al.l [20051: iTacconi et all I2006L 120081 : 
iNaravanan et al.l 12009 ). If the star formation continues 
with SFR ^ a few 10^ yr~^, the gas consumption 
time is ~ a few 10 Myr. This scenario is well described 
by hydrodynamic simulations of Naravanan et al. (2009 ) 
in which a major merger with a ^10^^ Mq dark matter 
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Fig. 5. — Synthetic SEDs of old stellar components (red), star- 
burst components at the age of 100 Myr with an initial SFR of 
3000 Mq yr"-*^ (blue), and a sum of both (green) created with the 
GALAXEV library. Photometry data are also shown. 



halo produces a peak flux of ~7-8 mJy at 850 /im for a 
short duration (<50 Myr). We are observing the short- 
lived glow of SXDF 850.6 as a bright SMG. Even if ah the 
molecular gas converts into stars, the total stellar mass 
will increase by at most 30%. This supports an idea that 
the large fraction of stellar mass in SMGs is constructed 
before the current submm bright phase (e.^.jBorvs et al.l 
[20051 : iDve etHI[2008l : 'Michalowski et al.""2009l 

The specific star formation rate (SSFR; SFR per unit 
stellar mass, an estimator of current star- forming activ- 
ity) of -5-20 Gyr-i places SXDF 850.6 on the high 
mass end of the correlation between stellar mass and SFR 
for z 2 star -forming galaxies or above the correlation 
(e.g lErb et a l. 2006b; Reddy et al. 2006; Daddi et al. 
l2007l ). If the large population of stars were constructed 
before the current starburst phase, the enhancement of 
star formation could be due to mergi ng of massive sta r 
forming galaxies (e.g., merging BzKsi iTakad et al.ll2008[ ). 



and it is possible that this galaxy evolves into a massive 
present-day elliptical through majo r merger events (e.g., 
iLillv et al-lllOOl iSmail et al.l[200l . 

7. SUMMARY 

We conducted SMA 880 /im observations of a submil- 
limeter galaxy, SXDF 850.6, detected with SCUBA in 
the Subaru/XMM-A/'ei6'ton Deep Field, which has mul- 
tiple radio and IR counterpart candidates. The precise 
astrometry with the SMA shows that the correct coun- 
terpart is the most distant radio source from the SCUBA 
coordinates. Although there is a lack of a corresponding 
source at the SMA position in deep images of 5, V ^ 
i\ z' bands, we find a counterpart at J, i^, 3.6 /im, 
4.5 /im, 5.8 /im, 8.0 /im, 24 /im, 1100 /im, and 21 cm. A 
detailed analysis of fitting a library of synthetic SEDs to 
the photometry from the ii-band to 8.0 /im gives the best- 
fit parameters of z = 1.87l^;J^, = 2.5toi x 10^^ 
Ay = 3.0t?;^ mag, and age = 720^2^0° Myr. The SED 
fit to the submm-radio photometry using the three dif- 
ferent starburst SED models provides the consistent red- 
shift of z - 1.8-2.5 and Lir = (7-26) x 10^^ Lq. The IR 
color diagnostic suggests that the AGN contribution to 
the SED is small. The SFR derived from the IR luminos- 
ity is 1300-4500 Mq yr-\ and the SSFR is -10 Gyr-^ 
If the intense starburst continues, the estimated molecu- 
lar gas mass of (2-5) x 10^^ Mq will last for a few Myr. 
Given these facts, it is suggested that SXDF 850.6 is a 
mature system that already contains a dominant frac- 
tion of its expected final stellar mass in an old stellar 
population, and which is also experiencing a significant 
starburst. The enhancement of star formation could be 
due to major mergers. 
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